Photon Filamentation in Resonant Media with High Presnel Numbers 
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The phenomenon of turbulent photon filamentation occurs in lasers and other active optical media 
at high Fresnel numbers. A description of this phenomenon is suggested. The solutions to evolution 
equations are presented in the form of a bunch of filaments chaotically distributed in space and having 
different radii. The probability distribution of patterns is defined characterizing the probabilistic 
weight of different filaments. The most probable filament radius and filament number are found, 
being in good agreement with experiment. 
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There exists a well established analogy between the evolution equations describing optical phenomena and hydrody- 
namic equations for compressible viscous liquid [1-3] . The Fresnel number for optical systems plays the same role as 
the Reynolds number for fluids. In optical systems, with increasing Fresnel number, there appears chaotic behaviour 
similar to fluid turbulence arising with increasing Reynolds number. Because of this analogy, one terms as the optical 
turbulence [1-3] the chaotic phenomena occurring in high-Fresnel-number optical systems. 

For an optical cavity of the standard cylindrical shape of radius R and length L, the Fresnel number is = nR^ /XL, 
with A being the optical wavelength. Physical processes accompanying the route from the regular behaviour to the 
turbulent regime, arising with increasing Fresnel numbers, are practically the same in different types of active optical 
media, such as lasers and photorefractive crystals. 

When Fresnel numbers are very small, F <C 1, there exists the sole transverse mode almost uniformly filling the 
optical cavity. At Fresnel numbers around F 1, several transverse modes emerge seen as bright spots in the 
transverse cross-section. These spatial structures arc regular both in space and in time, since in space they form 
ordered geometric arrays, as polygons, and because in time they are cither stationary or oscillating periodically. 
These transverse structures are imposed by the cavity geometry and correspond to the empty-cavity Gauss-Laguerre 
modes. The appearance of such regular geometric structures is well understood theoretically, with the theory being 
in good agreement with experiments both for lasers [4-9] as well as for photorefractive crystals [10-12]. For Fresnel 
numbers up to F w 5, the number of regular modes is proportional to F^. 

The transition from the regular behaviour to a turbulent state occurs around the Fresnel number F « 10. For 
Fresnel numbers F > 15, the arising spatial structures are principally different from those defined by the empty-cavity 
Gauss-Lagucrrc modes. The related modal description becomes no longer relevant and the boundary conditions do 
not play the role any more. In the medium, there appears a large number of filaments stretched along the cavity 
axis. In the transverse cross-section, they are distributed chaotically. The correlation length is much shorter than the 
mean distance between filaments, so that they radiate independently of each other, aperiodically flashing in time. The 
number of filaments, contrary to the regular regime, is proportional to F. This chaotic filamcntation has been observed 
in both photorefractive crystals [10-12], e.g. in Bii2Si02o, as well as in many lasers, such as Ne, Tl, Pb, N2, and 
vapor lasers [13-17], CO2 and dye lasers [18-26]. The specific random properties of the observed filamentation, such 
as chaotic distribution of filaments in space, the absence of correlations between them, and their aperiodic flashing, 
are analogous to the properties of turbulent fluids [1-3]. Because of this analogy of the chaotic photon filamentation in 
high-Fresnel-number optical media with turbulent phenomena in high-Reynolds-number viscous fiuids, one commonly 
calls [18-26] the former the optical turbulence or, because of the formation of bright filaments with a high density of 
photons, it can be called the turbulent photon filamentation. 

Contrary to the low- Fresnel number case, with regular optical structures, for which a good theoretical understanding 
has been achieved [4-12], for the high-Fresnel-number regime of the turbulent photon filamentation, there has been 
no persuasive theory suggested. This problem was addressed in Refs. [27-30]. However, because of the complexity of 
the phenomenon, only some oversimplified stationary models were considered. In the present paper, the description 
of the phenomenon of turbulent photon filamcntation is suggested, based on realistic evolution equations. It has 
become possible to solve now this problem owing to the recently developed method for treating nonlinear differential 
equations, called the scale separation approach [31-33]. 

For the system of N resonant atoms interacting with electromagnetic field, we use the standard dipole approxima- 
tion. The evolution equations to be considered are those for the averages of atomic operators, 

u{r,t) = <a-{r,t)>, s{f,t) = < a%r,t) > , (1) 

where a~ is a transition operator and is the population-difference operator. Equations for the quantities (1) can be 
easily obtained following the usual way, by writing the Heisenberg equations of motion, eliminating the field variables, 
and employing the standard semiclassical and Born approximations. All this machinery is well known and can be 
found, for instance, in the book [34]. In order to present resulting equations in a compact form, it is convenient 
to introduce the notation f = fo + f^ad for an effective field acting on an atom. This field consists of the term 
/o = —id ■ Eq, due to the cavity seed field, and of the term 

frad{r,t)= - \i-ip j [^p{r-r')u{r\t)-e-d^ ip*{r-r')u*{r',t)\ dr' (2) 

describing the effective interaction of atoms through the common radiation field. Here p is the density of atoms, 
d = doSd is the transition dipole, and ip{r) = exp(iA;o|f|)/A;o|f|, with fco = ojo/c and 7 = Ak^d^/S. The cavity seed 
field selecting a longitudinal propagating mode, is taken, as usual, in the form of plane waves, ^0 = ^le'^'^^"'^*^ -|- 
g*^-t{kz-wt) _ rpjjg resulting equations are 

fill H s 
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dt 



272H' + s(«*/ + r«) , (3) 



where loq is the transition frequency, 71 and 72 are the longitudinal and transverse relaxation widths, respectively, 
and C is the pumping parameter. 

We consider a cylindrical cavity, typical of lasers, with radius R and length L. For the case of high Fresnel numbers 
F » 10, the cavity can house several transverse modes. Therefore, it is reasonable to look for solutions of Eqs. (3) 
in the form of a bunch of Nf filaments, so that 

Nf Nf 

u{r,t)= ^ Un{f,t)Qn{x,y)e'^'', s{f,t) = ^ Sn{f,t)&n{x,y) , (4) 

n— 1 n— 1 

where 6„(x, y) = 6(6^ — (x — x„)^ + (?/ — is a unit-step function; the pair x„ and j/„ defines the axis of a filament; 
and bn is the radial distance from the filament axis, at which the solutions decrease by an order of magnit ude. If the 
filament profile is approximately of normal law exp(— r^/2r^J, with r„ being the filament radius, then 6„ = \/2 In 10 r„. 
The filaments do not interact with each other, because of which they cannot form a regular space structure, and the 
locations of their axes in the transverse cross-section are random. The absence of interactions between filaments is 
due to the fact that the kernel ip{r), playing the role of a Green function in the effective interaction field (2), fastly 
oscillates and diminishes with increasing r. Let us stress that the presentation of solutions in the form of a bunch 
of filaments (4) does not exclude the possible case of a sole filament almost uniformly filling the whole sample. This 
is because the number of filaments and their radii will be defined in a self-consistent way. Substituting the form (4) 
into Eqs. (3), we obtain a system of equations for «„, s„, and To simplify these equations, we employ the 

mean-field approximation for the averages 

u{t)= / Un{f,t)df, s{t)= / Sn{r,t)df, (5) 

''ra JVn JVn 

and for the similarly defined where the integration is over the volume F„ = nb'^L of a cylinder enveloping a 

filament. For the simplicity of notation, the index n labelling filaments is dropped from the left-hand sides of Eqs. 
(5). In this way, Eqs. (3) are transformed to the system of equations 

— = + T)u - isd ■ El e"''"* , 

J = -4gj2\u\^ - 71 (s - C) - 4 Im {u*d- Ei e"^-*) , (6) 



d\ul_ 
dt 



= -2T\uf + 2slm( u*d ■ Ei e'^* 



for the functions (5). Here the notations for the collective frequency ft = lvq + g'j2S and for the collective width 
r = 72(1 — gs) are used, with the effective coupling parameter 



472 K Jv^ ko\r 

and with g' defined as is Eq. (7) except for sine being replaced by cosine. 

Equations (6) describe the evolution of a filament. These evolution equations can be analysed by using the scale 
separation approach [31-33]. For applying the latter, we, first, take into account the existence of usual small parameters 
71 /wo -C 1 and 72 /wo 1, and also we consider the quasiresonance case, when |A|/a;o 1, with the detuning 
A = w — ujQ. Then the functional variable u in Eqs. (6) is to be classified as fast, as compared to the slow variables s 
and The latter are quasi-invariants for the fast variable. Following the scale separation approach [31-33], we solve 
the equation for the fast variable u, which is straightforward when s is a quasi-invariant. The found solution for u 
has to be substituted into the right-hand sides of the equations for the slow variables s and with averaging these 
right-hand sides over explicitly entering time of fast oscillations. This procedure, complimented by the transformation 



3 



where a is assumed to be small, a <C 1, leads to the equations 



ds dw 

— = -4572a' - 71 (s - C) , — = -2-f2{l - gs)w (9) 

for the slow variables. 

The properties of solutions to Eqs. (9) essentially depend on the value of the effective coupling g defined in Eq. 
(7). From this definition, it is evident that the value of g depends on the characteristics of the considered filament, 
in particular, on the filament radius r„ which enters Eq. (7) through the radius 6„ of the enveloping cylinder. It 
is convenient to introduce the dimensionless parameter (3 = nb^/XL and to consider g = g{/3) as a function of this 
parameter in the domain < /? < F. Thus, each filament can be characterized by its radius or, equivalently, by the 
related parameter p. The bunch of filaments parametrized by /3 forms the solutions (4). The stability properties of 
filaments with different values of fi are different, which can be described by defining the probability distribution of 
patterns corresponding to filaments with /3 varying in the interval [0,F]. The latter can be done as follows. Any 
system of evolution equations, like Eqs. (9), can always be presented in the normal form dy/dt = v, in which y = {yi}, 
and V = {vi} is a velocity field. In the case of Eqs. (9), we have yi = s, y2 = w. 

As is known from statistical mechanics, a probability p is connected with entropy S by the relation p ~ . In 
the nonequilibrium case, it is more appropriate to count the entropy S{t), which is a function of time, from its initial 
value 5(0), thus, considering the entropy variation AS{t) = S{t) — 5(0). The entropy is defined as the logarithm of an 
elementary phase volume, S{t) = ln|A$(t)|, where, in the nonequilibrium case, A$(t) = YljSyj{t). The elementary 
phase volume can be expressed as A$(f) = ]J - J2j -^ii(*)<52/j(0) through the elements Mij{t) = 5yi{t)/Syj{0) of the 
multiplier matrix M{t) = [Mjj(t)]. Hence, the entropy variation is A5(t) = TrL, with the matrix L = [Ljj] being 
composed of the elements Lij = InM^ . By variating the evolution equations, one gets the equation dM/dt = JM 
for the multiplier matrix M{t), where J = [Jij] is the Jacobian matrix with the elements Jij = Svi/Syj. In this way, 
the probability p ~ e^"^ acquires the form p ~ exp{— Tr£}. From the equation for the multiplier matrix, one has 
dTiL/dt — TrJ, so that the entropy variation is AS{t) = A{t')dt' , where A = Tr J is called the contraction rate 
[35]. Therefore, for the probability distribution of patterns, labelled by a parameter /3, one has 
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p{(3,t)= ^ expj- A{f3,t')dt'j , (10) 

where Z is the normalization factor. As is clear, that pattern is preferred over the other which possesses a higher 
probability weight. This implies, because of the form (10), that the most preferable pattern is that for which the 
contraction is minimal. At the initial stage of pattern formation, one may write p ^ exp{— A(/3, 0)t}. Hence, the most 
probable pattern to be formed is that corresponding to the minimal contraction rate A(/3,0). 

Considering the evolution equations (9), we keep in mind the standard way of exciting laser media by means of a 
nonresonant pumping described by the pumping parameter C > 0. No resonant fields are involved, so that s(0) < 0. 
It is easy to calculate the contraction rate A(/3, 0) corresponding to Eqs. (9), which is A(/3,0) = —71 — 272(1 — gso), 
where sq = s(0). The minimum of this contraction rate is provided by the maximum of g = g{P). The maximum of 
g{P), with g{(3) given by Eq. (7), occurs at /3 = 0.96. From here we find the most probable filament radius 



r/ = 0.26VAL . 

The most probable number of filaments can be evaluated from the normalization integral y J s{f,t)df = where 
the integration runs over the whole volume of the sample, V = ttR^L. Considering this normalization integral for the 
stage when the filaments have already been formed and the population difference inside a filament of radius r/ has 
reached a value close to (, we obtain 



The number of filaments is proportional to the Presnel number F, which is common for the turbulent photon filamen- 
tation. 

Calculating the filament radius r/ for the variety of experiments with different lasers, we find that the theoretical 
values of r/ are in perfect agreement with all experimental data available. Thus, for Ne, Tl, Pb, N2, and N2' laser 
[13-17], we find rj w 0.01 cm and for CO2 and dye lasers [18 24], we have rj w 0.08 cm and rj w 0.01 cm, respectively. 

In conclusion, the solution for the problem of turbulent photon filamentation occurring in resonant media at high 
Presnel numbers is suggested. The consideration is based on realistic evolution equations for resonant atoms under 
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conditions typical of lasers. The solutions to these equations have the form corresponding to a bunch of filaments with 
different radii. The probability distribution of filament radii is found self-consistently from the evolution equations. 
The results are in good agreement with all experiments on the turbulent photon filamentation in laser media [13 24]. 

This work has been supported by the Bogolubov-Infeld Grant of the State Agency for Atomic Energy, Poland and 
by a Grant of the Sao Paulo State Research Foundation, Brazil. 
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